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ABSTRACT 

A big challenge in solar and stellar physics in the coming years will be to de- 
cipher the magnetism of the solar outer atmosphere (chromosphere and corona) 
along with its dynamic coupling with the magnetic fields of the underlying pho- 
tosphere. To this end, it is important to develop rigorous diagnostic tools for 
the physical interpretation of spectropolarimetric observations in suitably chosen 
spectral lines. Here we present a computer program for the synthesis and inver- 
sion of Stokes profiles caused by the joint action of atomic level polarization and 
the Hanle and Zeeman effects in some spectral lines of diagnostic interest, such 
as those of the He I 10830 A and 5876 A (or D3) multiplets. It is based on the 
quantum theory of spectral line polarization, which takes into account in a rigor- 
ous way all the relevant physical mechanisms and ingredients (optical pumping, 
atomic level polarization, level crossings and repulsions, Zeeman, Paschen-Back 
and Hanle effects). The infiuence of radiative transfer on the emergent spectral 
line radiation is taken into account through a suitable slab model. The user can 
either calculate the emergent intensity and polarization for any given magnetic 
field vector or infer the dynamical and magnetic properties from the observed 
Stokes profiles via an efficient inversion algorithm based on global optimization 
methods. The reliability of the forward modeling and inversion code presented 
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here is demonstrated through several apphcations, which range from the infer- 
ence of the magnetic field vector in solar active regions to determining whether 
or not it is canopy-like in quiet chromospheric regions. This user-friendly diag- 
nostic tool called "HAZEL" (from HAnle and ZEeman Light) is offered to the 
astrophysical community, with the hope that it will facilitate new advances in 
solar and stellar physics. 

Subject headings: magnetic fields — polarization — radiative transfer — scatter- 
ing — Sun: chromosphere — methods: data analysis, numerical 



Introduction 



The present paper describes a computer program for the synthesis and inversion of 
Stokes profiles resulting from the joint action of the Hanle and Zeeman effects in some 
spectral lines of diagnostic interest, such as those pertaining to the He 1 10830 A and 5876 A 
(or D3) multiplets. The effects of radiative transfer on the emergent spectral line radiation are 
taken into account through a suitable slab model. Our aim is to provide the solar and stellar 
physics communities with a robust but user-friendly tool for understanding and interpreting 
spectropolarimetric observations, with the hope that this will facilitate new advances in solar 
and stellar physics. 

In particular, the lines of neutral helium at 10830 A are of great interest for empirical 
investigations of the dynamic and magnetic proper ties of plasma structures in the solar 



chromosphere arid corona, such as acti ve regions fe.g.. lHarvev fc 1 



Lagg et al. 2004; Centeno et al. 20061), filamen ts (e.g., Lin et al] 



2002bh. 



prommences (e.g. , 



all 



1971 



19981: 



Merenda et al.ll2006l) and spicules (e.g.. lTrujillo Bueno et al.ll2005l : 



Riiedi et al 



1996 



Truiillo Bueno et al. 



Socas-Navarro fc Elmorell2005l ). The same applies to the lines of the He I D3 multiplet at 
5876 A which have been used for investigating the magnetic field vect o r in solar prominences 



and spicules (e.g. , iLandi DegFInnocenti 



Casini et al.ll2003l : iLopez Ariste fc Casini 



1982 



2005 



Querfeld et a 


1985 


■ 


Bommier et al. 


Ramelli et al. 


2006a 


b 


) 



Such helium lines result from transitions between terms of the triplet system of helium 
(ortho-helium), whose respective J-levels (with J the level's total angular momentum) are 
far less populated than the ground level of helium (that is, than the singlet level ^Sq), except 
perhaps in flaring regions. On the other hand, the lower term (2s^Si) of the He I 10830 A 
multiplet is the ground level of ortho-helium, while its upper term (2p'^P2,i,o) is the lower one 
of 5876 A (whose upper term is 3d'^D3 2,1). Therefore, the significant difference in the ensuing 
optical thicknesses of the observed solar plasma structure implies that when the radiation 
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in these spectral lines is observed on the solar disk it is much easier to see structures in 
10830 A than in 5876 A, while both lines are clearly seen in emission when observing off- 
limb structures such as prominences and spicules. The additional fact that the Hanle effect 
in forward scattering creates measurable linear polarization signals in the lines of the He I 
10830 A multiplet when the magnetic field is inclined with respect to the local solar vertical 



direction (jTrujillo Bueno et al.l l2002bl ) . and that there is a nearby photospheric line of Si I, 
makes the 10830 A spectral region very suitable for investigating the coupling between the 
photosphere and the corona. 

While the Stokes / profiles of the 10830 A and 5876 A helium lines depend mainly on the 
distribution of the populations of their respective upper ( J^) and lower (J/) levels along the 
line-of-sight (LOS), their Stokes Q, U and V profiles depend on the strengths and wavelength 
positions of the vr (AM = M„ -M/ = 0), (Xbiue (AM = +1) and a^ed (AM = -1) transitions, 
which can only be calculated correctly within the framework of the Paschen-Back effect 
theory. Moreover, the Q, U and V profiles are also affected by the atomic level polariz ation 



induced by anisotropic pumping processes (e.g.. lLandi DeglTnnocenti fc Landolfill2004j ). An 
atomic level of total angular momentum J is said to be polarized when its magnetic sublevels 
are unequally populated and/or when there are quantum coherences between them. The 
radiative transitions induced by the anisotropic illumination of the helium atoms in the 
solar atmosphere are able to create a significant amount of atomic pola rization in the helium 



levels, even in the metastable lower level of the He I 10830 A multiplet ( jTrujillo Bueno et al. 



2002bl ). If the net circular polarization of the incident radiation at the wavelengths of the 
helium transitions is negligible, as it uses to be the case, the radiatively induced atomic 
level polarization is such that the populations of substates with different values of \M\ 
are different (non-zero atomic alignment), while substates with magnetic quantum numbers 
M and — M are equally populated (zero atomic orientation). On the other hand, elastic 
collisions with the neutral hydrogen atoms of the solar chromospheric and coronal structures 
are unable to destroy the atomic polarization of the He I levels. As a result, even in the 
absence of magnetic fields, linearly polarized spectral line radiation would be produced, 
simply because the population imbalances among the magnetic sublevels imply more or 
fewer vr-transitions, per unit volume and time, than a transitions. The atomic polarization 
of the upper level of the line transition under consideration is thus responsible of a selective 
emission of polarization components, while that of the lower level may give rise to a selective 
absorption of polarization components ( "zero- field" dichroism). In order for this type of 
dichroism to produce a measurable contribution to the emergent linear polarization it is 
necessary to have a substantial line-center optical thickness along the LOS or, assuming a 
small but non-negligible optical thickness, that the plasma st ructure under consideration 



is observed against the bright background of the solar disk (jTrujillo Bueno et al.l l2002bl : 
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Trujillo Bueno &: Asensio Ramoa 120071 ) . Therefore, the observable effects of dichroism are 
easier to detect in He I 10830 A than in 5876 A. 

In the presence of a magnetic field the emergent polarization changes because of the 
following two reasons. First, because a magnetic field modifies the atomic level polarization, 
not only by producing the Hanle-effect relaxation of the quantum coherences pertaining to 
each individual J-level, but also through possible interferences between the magnetic sub- 
levels pertaining to different J-levels, which give rise to a variety of remarkable effects such 
as the transfer of atomic a lignment to atomic orientat ion in the J-levels of the upper term 



of the He I D3 multiplet (iLandi Degl'Innocentil Il982l ) or the enhance ment of the scatter- 



ing po larization in the D2 line of Na I by a vertical magnetic field (iTrujillo Bueno et al. 



2002al ). Second, because the magnetic splitting of the atomic energy levels give rise to signif- 
icant wavelength shifts between the vr and a transitions (as compared with the spectral line 
width) and, consequently, to the generation of measurable linear and/or circular polariza- 
tion (i.e., the familiar transverse and longitudinal Zeeman effects, respectively). Obviously, 
a correct modeling of the spectral line polarization that results from the joint action of the 
Hanle and Zeeman effects requires the application of the quantum theor y of spectral line 
polarization (see the monograph by lLandi Degl'Innocenti fc Landolf]||2004j ). as done by sev- 
eral researchers for interpreti ng spectropolarimetric observations of solar plasma structures 



2OO3I: iLopez Ariste fc Casin 



2002bl . 



2005 



Merenda et al. 



in the He I D,s multiplet fe.g..lLandi Degl'Innocentilll982l : iBommier et al 



20061 ) 



2005) and in the He 1 10830 A triplet (e.g.. ITrujillo Bueno et al. 



1994 



Casini et al. 



Over the last few years new computer programs for the synthesis and inversion of Stokes 
profiles induced by the joint action of atomic level polarization and the Hanle and Zee- 
man effects have been devel oped and applied to the interpre t ation of spectropolarimetric 



observations. For example, iLandi DegPInnocenti fc Landolfil (120041 ) have developed some 



forward modeling codes with which they have calculated several Hanle effect diagrams and 
theoretic al Stokes profiles of lines fro m complex atomic models, while the computer pro- 



gram of ICasini fc Manso Saina (120051 ) can treat even the case of a hyperfine structured 



multiplet taking into account the quantum interferences between the F levels belonging to 



that 


Casini et al. 


( 


2003 


) api 


(Lopez Ariste & Casini 


2OO2I 



to spectropolarimetric observations of solar prominences in the 
He I D3 multiplet, providing two-dimensional maps of the magnetic field vec tor and further 



evidence for strengths significantly larger than average. On the other hand, iMerenda et al. 



( 20061 ) opted for an inversion strategy for the He I 10830 A multiplet in which the longitudi- 
nal component of the magnetic field vector is obtained from the measured Stokes V profiles 
(which are dominated by the Zeeman effect for the lines of the He I 10830 A multiplet), and 
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its orientation from the observed Stokes Q and U profiles (which in solar prominences are 
due to the presence of atomic level polarization). In the just mentioned computer programs 
and applications the optically thin assumption was used, which was a suitable approximation 
for the particular prominences observed, but an unreliable one in general (and especially for 
the interpretation of observations in the He I 10830 A multiplet). 

A simple but suitable model for taking into account r adiative tran sfer effects is the 
constant-property slab model used by iTrujillo Bueno et al.l (l2002bl . l2005l ) for the interpre- 
tation of spectropolar imetric observations in solar filaments a nd spicules, which has been 
recently extended by ITrujillo Bueno &: Asensio Ramod (120071 ) to include magneto-optical 
effects. These authors have applied this "cloud" model for the interpretation of spectropo- 
larimetric observations in order to point out that the atomic polarization of the helium 
levels may have an important impact on the emergent linear polarization of the He I 10830 
A multiplet, even for magnetic field strengths as large as 1000 G. Therefore, inversion codes 
th at neglect the infl uenc e of atomic level pol a rizati on, such as the Milne-Eddington codes 
of iLagg et al.l (120041 ) and ISocas-Navarro et al.l (120041 ). should ideally be used only for the in- 
version of Stokes proflles emerging from strongly magnetized regions (with B > 1000 G), or 
when the observed Stokes Q and U proflles turn out to be dominated by the transverse Zee- 
man effect (e.g., as it happens with s ome active regions fllaments as a result of the particular 
illumination conditions explained in iTrujillo Bueno &: Asensio Ramosll2007l ). 



It is also necessary to mention that iManso Sainz &: Trujillo Buend (l2003al ) developed 
a general radiative transfer computer program for solving the so-called non-LTE problem 
of the 2nd kind -that is, the multilevel scattering polarization problem including the Hanle 
effect of a weak magnetic held. These authors considered the multilevel atom approximation 
(which neglects quantum interferences between the sublevels of different J-levels), but took 
fully into account the effects of radiative transfer in realistic atmospheric models and the role 
of elastic and inelastic collisions in addition to all the relevant optical pumping mechanisms. 
An inter esting application using a sera i -empiri cal model of the solar atmosphere can be 
found in iManso Sainz &: Trujillo Buend (l2003bl ). The recent advances in the development 
of efficient numerica l methods for the so lution of non-LTE polarization transfer problems 
(e.g., the review by iTrujillo Buend 120031 ) and in computer technology make now possible 
even the performing of three-di mensional radiative transfer simul ations of the Hanle effect 
in convective atmospheres (e.g.. iTrujillo Bueno fc Shchukinal 120071 ). 



The previous introductory paragraphs strongly suggest that it would be of great interest 
to develop a robust but user-friendly computer program for the synthesis and inversion of 
Stokes profiles resulting from the joint action of atomic level polarization and the Hanle 
and Zeeman effects. We have done this by implementing an efficient global optimization 
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method for the solution of the inversion problem and by calculating at each iterative step 
the emergent spectral line polarization through the solution of the Stokes-vector transfer 
equation in a slab of constant physical properties in which the radiatively-induced atomic 
level polarization is assumed to be dominated by the photospheric continuum radiation. At 
each point of the observed field of view the slab's optical thickness is chosen to fit the Stokes 
/ profile, which is a strategy that accounts implicitly for the true physical mechanisms that 
populate the triplet levels of He I (e.g., the photoionization-recombination mechanism, as 



shown by lCenteno et al.ll2008l ). The observed Stokes Q, U and V profiles are then used to 



infer the magnetic field vector along with a few extra physical quantities. 

The outline of this paper is as follows. The formulation of the problem is presented 
in §2, where we review the relevant equations within the framework of the quantum theory 
of spectral line polarization. The forward modelling option of our computer program is 
described in §3, including some interesting examples of possible applications. §4 deals with a 
detailed description of the Stokes inversion option, while §5 considers a variety of applications 
aimed at an in-depth testing of this diagnostic tool. The important issue of the possible 
ambiguities and degeneracies is addressed in §6, with emphasis on the Van Vleck ambiguity 
and on the possibility of inferring the atmospheric height at which the observed on-disk 
plasma structure is located. Finally, in §7 we summarize the main conclusions and comment 
on some ongoing developments. 



2. Formulation of the problem 

We consider a constant-property slab of He I atoms, located at a height h above the 
visible solar "surface" , in the presence of a deterministic magnetic field of arbitrary strength 
B, inclination 9b and azimuth xb (see Fig. 1). The slab's optical thickness at the wavelength 
and line of sight under consideration is r. We assume that all the He I atoms inside this slab 
are illuminated from below by the photos pheric so lar continuum radiation field, whose center- 



to-limb variation has been tabulated by iPiercd (|2000| ). The ensuing anisotropic radiation 



pumping produces population imbalances and quantum coherences between pairs of magnetic 
sublevels, even among those pertaining to the different J-levels of the adopted He I atomic 
model. This atomic level polarization and the Zeeman-induced wavelength shifts between 
the 71 (AM = Mu- Ml = 0), cxbiue (AM = +1) and a^ed (AM = -1) transitions produce 
polarization in the emergent spectral line radiation. 

In order to facilitate the reading and understanding of this paper, in the following 
subsections we summarize the basic equations which allow us to calculate the spectral line 
polarization taking rigorously into account the joint action of the Hanle and Zeeman effects. 
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To this end, we have apphed the quantum t heory of spectral hne polarization, wh ich is 
described in great detail in the monograph by lLandi DegFInnocenti fc Landolfil (120041 ) . We 
have also applied several methods of solution of the Stokes-vector transfer equation, some of 
which can be conside red as particular cases of the two general methods explained in §6 of 



Truiillo Buenol (120031 ). 



2.1. The radiative transfer approach 

The emergent Stokes vector I(z/, f2) = (J, Q, U, (with t=transpose, z/ the frequency 
and fl the unit vector indicating the direction of propagation of the ray) is obtained by 
solving the radiative transfer equation 



as 



e(z/, n) - K(z/, n)i{u, n), 



where s is the geometrical distance along the ray under consideration, e(i/, fl) = (e/, eg, eu, ey)^ 
is the emission vector and 
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Vi / 



(2) 



is the propagation matrix. Alternatively, introducing the optical distance along the ray, 
dr = —rjids, one can write the Stokes-vector transfer Eq. ([1]) in the following two ways: 



The first one, whose formal solution requires the use of the ev olution operator intro- 
duced by lLandi Degl'Innocenti fc Landi Degl'Innocenti I (119851 ). is 



K*l - S, 



(3) 



where K* = K/r?;- and S = e/rir . The formal solution of this equation can be seen in 
eq. (23) of iTruiillo Buenol J2003h . 



The second one, whose forr nal solution does n ot require the use of the above-mentioned 
evolution operator is (e.g.. lRees et al. 1119891 ) 



dr 



5eff, 



(4) 
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where the effective source-function vector Sefr = S — K'I, being K' = K* — 1 (with 
1 the unit matrix). The formal solution of this equation can be seen in eq. (26) of 



Truiillo Buenol fl2003f ). 



Once the coefficients e/ and ex (with X = Q,U,V) of the emission vector and the 
coefficients rji, rjx, and px of the 4x4 propagation matrix are known at each point within 
the medium it is possible to solve formally Eq. ([3]) or Eq. (jl]) for obtaining the emergent 
Stokes profiles for any desired line of sight. Our computer program considers the following 
levels of sophistication for the solution of the radiative transfer equation: 



Numerical Solutions. The most general case, where the properties of the slab vary 
along the ray path, has to be solved numerically. To this end, two efficient and accu- 
rate methods o f solution of the Stokes- vector transfer equation are those proposed by 
Trujillo Buenol (120031 ) (see his eqs. (24) and (27), respectively). The starting points 
for the development of these two numerical methods were Eq. ([3]) and Eq. (jl]), re- 
spectively. Both methods can be considered as generalizations, to the Stokes-vector 
transfer case, of the well-known short characteristics method for the solution of the 
standard (scalar) transfer equation. 

Exact analytical solution of the problem of a constant-property slab including the magneto- 
optical terms of the propagation matrix. For the general case of a constant-property 



slab of arbitrary optical thickness we actually have the followi ng analytical s o 



which can be easily obtained as a particular case of eq. (24) of iTrujillo Buenol (120031 ): 



ution. 



I 



e-^*"I 



(5) 



where Isun is the Stokes vector that illuminates the slab's boundary that is most distant 
from the observer. We point out that the exponenti al of the propagation matrix K* 
has an analytical expression similar to eq. (8.23) in iLandi Degl'Innocenti fc Landolfi 



Approximate analytical solution of the problem of a constant-property slab including the 
magneto-optical terms of the propagation matrix. An approximate analytical solution 
to the con stant-property slab pr oblem can be easily obtained as a particular case of 
eq. (27) of iTruiillo Buenol J2003h : 



I = [1 + ^oK'] ' [{e-^l - ^mK') Is^n + (^M + ^o)S] , 



(6) 
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where the coefficients '^m and \l/o depend only on the optical thickness of the slab at 
the frequency and line-of-sight under consideration, since their expressions are: 



^0 



1 - 



— e 



T 



(7) 



Note that Eq. for the emergent Stokes vector is the one used by lTrujillo Bueno fc Asensio Ramos 
(120071 ) for investigating the impact of atomic level polarization on the Stokes profiles 
of the He I 10830 A multiplet. We point out that, strictly speaking, it can be con- 
sidered only as the exact analytical solution of the optically-thin constant-property 
slab problenj^. The reason why Eq. ([6]) is, in general, an approximate expression 
for calculating the emergent Stokes vector is because its derivation assumes that the 
Stokes vector within the slab varies linearly with the optical distance. However, it 
provides a fairly good approximation to the emergent Stokes profiles (at least for all 
the problems we have investigated in t his paper). Moreover, the results of fig. 2 of 
Trujillo Bueno fc Asensio Ramos! (120071 ) remain also virtually the same when using in- 
stead the exact Eq. ([5]), which from a computational viewpoint is significantly less 
efficient than the approximate Eq. ([6]). 

Exact analytical solution of the problem of a constant-property slab when neglect- 
ing the second-order terms of the Stokes-vector transfer equation. Simplified expres- 
sions for the emergent Stokes vector can be obtained when e/^ex and rji^irjx, Px)i 
which justifies to neglect the second-order terms of Eq. ([1]). The resulting approxi- 
mate formulae for the emergent Stokes parameters are given by eqs. (9) a nd (10) of 



Truiillo Bueno &: Asensio Ramod (120071 ). which are identical to those used by lTrujillo Bueno et al, 
(l2005h for modeling the Stokes profiles observed in solar chromospheric spicules. We 
point out that there is a typing error in the sen tence that introduces such eqs. (9) and 
(10) in lTrujillo Bueno fc Asensio Ramos! (12007! ). since they are obtained only when the 
above-mentioned second-order terms are neglected in Eq. ([1]), although it is true that 
there are no magneto-optical terms in the resulting equations. 

Optically thin limit. Finally, the most simple solution is obtained when taking the 
optically thin limit (r^l) in th e equations reported in the previous poin t, which lead 
to the equations (11) and (12) of !Trujillo Bueno fc Asensio Ramos! (12007! ). Note that if 

Isii 



(i.e., Jo = Xq = 0), then such optically thin equations imply that Xjl^txl^i- 



^More precisely, when the optical thickness of the slab is small in comparison with the eigenvalues of the 
matrix K'. 
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The coefficients of the emission vector and of the propagation matrix depend on the 
multipolar components, Pq{J,J'), of the atomic density matrix. Let us recall now the 
meaning of these physical quantities and how to calculate them in the presence of an arbitrary 
magnetic field under given illumination conditions. 



2.2. The multipolar components of the atomic density matrix 

We quantify the atomic polarization of the He I levels using the multipolar components 
of the atomic density matr ix. The He I atom can be correctly described under the framework 



of the L-S coupling (e.g., ICondon &: Shortleylll935l ). As illustrated in Fig. 2, the different 



J-levels are grouped in terms with well defined values of the electronic angular momentum 
L and the spin S. Since the '^He atoms are devoid of nuclear angular momentum, we do not 
have to consider hyperfine structur^. The energy separation between the J-levels pertaining 
to each term is very small in comparison with the energy difference between different terms. 
For example, the energy separation between the J = 3 and J = 2 levels of the term 3d^D (the 
upper term of the D3 multiplet) is of the order of 0.0025 cm~^, which is ~2xl0^ times smaller 
than the separation between the 3d^D and 3p^P terms. On the other hand, although the 
energy separations between the J-levels of the upper terms of the 10830 A and D3 multiplets 
are much larger than their natural widths, such J-levels suffer crossings and repulsions for 
the typical magnetic strengths encountered in the solar atmospheric plasma (e.g., the J = 2 
and J = 1 levels of the upper term of the He 1 10830 A multiplet cross for magnetic strengths 
between 400 G and 1600 G, while the J = 3 and J = 2 levels of the upper term of the D3 
multiplet cross for strengths between 10 G and 100 G, approximately). This can be seen 
clearly in Fig. [3l Therefore, it turns out to be fundamental to allow for coherences between 
different J-levels pertaining to the same term but not between the J-levels pertaining to 
different terms. As a re sult, we can represent the atom under the formalism of the multi- 



term atom discussed by lLandi Degl'Innocenti &: Landolfil (120041 ). 



In the absence of magnetic fields the energy eigenvectors can be written using Dirac's 
notation as \PLSJM), where /? indicates a set of inner quantum numbers specifying the 
electronic configuration. In general, if a magnetic field of arbitrary strength is present, the 
vectors \(3LSJM) are no longer eigenf unctions of the total Hamiltonian and J is no longer 
a good quantum number. In this case, the eigenfunctions of the full Hamiltonian can be 



^See iBelluzzi. Truiillo Bueno. &: Landi Degl'Innocentil (|2007l ) for the formulatfon and solution of an in- 
teresting problem where hyperfine structure is important. 
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written as the following linear combination: 

\f3LSjM) = J2c^ji(^LS,M)\f3LSJM), (8) 
J 

where j is a pseudo-quantum number which is used for labeling the energy eigenstates 
belonging to the subspace corresponding to assigned values of the quantum numbers P, 
L, S, and M, and where the coefficients Cj can be chosen to be real. 

In the presence of a magnetic field sufficiently weak so that the magnetic energy is 
much smaller than the energy intervals between the J-levels, the energy eigenvectors are 
still of the form \(3LSJM) {Cj{PLS, M) ^ and the splitting of the magnetic sublevels 
pertaining to each J-level is linear with the magnetic field strength. For stronger magnetic 
fields, we enter the incomplete Paschen-Back effect regime in which the energy eigenvectors 
are of the general form given by Eq. ([8]), and the splitting among the various M-sublevels is 
no longer linear with the magnetic strength. This regime is reached for magnetic strengths 
of the order of 10 G for the He I D3 multiplet and of the order of 100 G for the 10830 A 
multiplet (see Fig. [3]). If the magnetic field strength is further increased we eventually reach 
the so-called complete Paschen-Back effect regime, where the energy eigenvectors are of the 
form \LSMlMs) and each L-S term splits into a number of components, each of which 
corresponding to particular values of {M^ + 2M5) . 

Within the framework of the multi-term atom model the atomic polarization of the 
energy levels is described with the aid of the density matrix elements 

p^^^(jM,j'M') = {^LSjMWLSj'M'), (9) 

where p is the atomic density matrix operator. Using the expression of the eigenfunctions of 
the total Hamiltonian given by Eq. ([S]), the density matrix elements can be rewritten as: 

p^^^{jM,fM') = CiiPLS, M)C^j,{(3LS, M')p^^^{JM, J'M'), (10) 
JJ' 

where p^^^{JM, J'M') are the density matrix elements on the basis of the eigenvectors 



\(3LSJM). 



Following iLandi Degl'Innocenti &: Landolfil (120041 ). it is helpful to use the spherical sta- 



tistical tensor representation, which is related to the previous one by the following linear 
combination: 



^^^p5(J,J') = C'jiPLS,M)Ci,iPLS,M') 

MM' 

-iy-'^V2KTll^l^ _^)p^^"(jM,/M'), (11) 



jj'MM' 

X 
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where the 3-j symbol is defined as indicated by any suitable textbook on Racah algebra. This 
alternative representation has some advantages. Firstly, the well-known results obtained 
when atomic polarization effects are disregarded are easily recovered by considering only 
the elements of the density matrix with K = and Q = 0. Secondly, the transformation 
law under rotations is much simpler because it involves only one rotation matrix. Finally, 
each Pq element has a clear physical interpretation: the pg elements (with Q = 0, ±1, ±2) 
are called the atomic alignment components, while the pg elements (with Q = 0, ±1) are 
the atomic orientation components. The ^^^Pq{J, J') elements are, in general, complex 
quantities (except for the elements with Q = and J = J', that are real quantities), so 
that, taking into account that the density matrix is an hermitian operator, the number of 
real quantities required to describe the atomic polarization properties of a given L-S term 
is (25* + 1)^(2L + 1)^. This makes a total of 405 real quantities to describe the atomic 
polarization of the 5-term model atom shown in Fig. [21 



2.3. Statistical equilibrium equations 



In order to obtain the Pq{J, J') elements we have to solve the statistical equilibrium 
equations. These equations, written in a reference system in which the quantization axis {Z) 
is directed along the magnetic field vector and ne glecting the infiuence of collisions, can be 
written as (iLandi DeglTnnocenti &: Landolfil 120041 ): 



|^^^p5(J,J') = -2'n\Y,Y.^^^s{KQJJ\K'Q'J"J"T'p^:{J\J"') 

K'Q' J" J'" 

+ ^ PeL.SpK.p^^ j'^)Ta{(3LSKQJJ', f3eLeSKeQeJeJ'e) 
+ ^^"^^""pQliJu, J'u) \fE{(3LSKQJJ\ f3^L^SKuQuJuJ'u) 

Lu QuJuJ'u 

+Ts{(3LSKQJJ', f3^L^SK^QuJuJ'J 
Yl ''''^Pq'{J'\J"')Wa{PLSKQJJ'K'Q'J"J"') 

K'Q' J" J'" 

+Re{PLSKQJJ'K'Q'J"J"') + Rs{(3LSKQJJ'K'Q'J"J"')] . (12) 

The first term in the right hand side of Eq. (fT2l) takes into account the infiuence of the 
magnetic field on the atomic level polarization. This ter m has its simplest expression in 
the c hosen magnetic field reference frame (see eq. 7.41 of iLandi Degl'Innocenti fc Landolfi 
20041 ). In any other reference system, a more complicated expression arises. The second. 
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third and fourth terms account, respectively, for coherence transfer due to absorption from 
lower levels (T^), spontaneous emission from upper levels (T^;) and stimulated emission 
from upper levels (T5). The remaining terms account for the relaxation of coherences due to 
absorption to upper levels (M^), spontaneous emission to lower levels (M^;) and stimulated 
emission to lower levels (M5), respectively. 

The stimulated emission and absorption transfer and relaxation rates depend explicitly 
on th e radiation field properties (see eqs. 7.45 and 7.46 of iLandi Degl'Innocenti fc Landolfi 
2OO4J ). The symmetry properties of the radiation field are accounted for by the spherical 
components of the radiation field tensor: 



3 

i=0 



(13) 



The quantities Tq^ii, ^) are spherical tensors that depend on the reference frame and on the 
ray direction O. They are given by 



(14) 



where R' is the rotation that carries the reference system defined by the line-of-sight f2 
and by the polarization unit vectors ei and 62 into the r eference system of the magnetic 
field, while Vpn(R') is the usual rot ation matrix (e.g., lEdmondj llQGOl ). Table 5.6 in 



Landi Degl'Innocenti fc Landolfi! (j2004j ) gives the Tq"{ 
eter Si (with Sq = I, Si = Q, S2 = U and S3 = V). 



O) values for each Stokes param- 



2.4. Emission and absorption coefficients 



Once the multipolar components ^^^Pq{J, J') are known, the coefficients e/ and ex (with 
X = Q,U,V) of the emission vector and the coefficients rjj, rjx, and px of the propagation 
matrix for a given transition between an upper t erm {l3L,i,S) and an lower term {0LpS ) 
can be calculated with the expressions of §7.6.b in lLandi Degl'Innocenti &: Landolfil (120041 ). 
These radiative transfer coefficients are proportional to the number density of He I atoms, 
jV. Their defining expressions contain also the V oigt profile and the Faraday- Voigt profile 
(see §5.4 in lLandi Degl'Innocenti fc Landolfil l2004l ). which involve the following parameters: 
a (i.e., the reduced damping constant), Wth (i.e., the velocity that characterizes the thermal 
motions, which broaden the line profiles), and fmac (i-e., the velocity of possible bulk motions 
in the plasma, which produce a Doppler shift). 
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It is important to emphasize that the expressions for the emission and absorption coef- 
ficients and those of the statistical equihbrium equations are written in the reference system 
whose quantization axis is parallel to the magnetic field. The following equation indicates 
how to obtain the density matrix elements in a new reference system: 

r'p'Q(J-J')]„ = E r'pi^-U'JXuKQcr, (15) 

Q' 

where Vq,q{R)* is the complex conjugate of the rotation matrix for the rotation R that 
carries the old reference system into the new one. 



2.5. The atomic model 

The atomic model we have used in our calculations includes the following five terms of 
the triplet system of neutral helium: 2s^S, 3s^S, 2p'^P, 3s^P and 3d^D (see Fig. [2]). It has been 
concluded that the inclusion of these five terms is sufficient for a reliable calculation of the 
atomic polarization that the anisotropic pumping of the ph otospheric cont inuum radiation 



produces in the lower and upper J-levels of the D3 multiplet (jBommierlll980l ) . Since the only 
level (with J = 1) of the lower term 2s^S is metastable, the adopted atomic model should 
be also satisfactory for the 10830 A multiplet. In order to check this we have compared the 
results obtained using two different atomic models. The first one is that of Fig. [2l The 
second one is a simplified version in which only the 2s^S, 2p^P and 3d'^D terms are taken 
into account. We have verified, for a large number of possible configurations, that differences 
in the resulting values of the multipolar components of the atomic density matrix are never 
larger than ~5% for the terms involved in the 10830 A transitions. Although the selection 
rules allow six transitions among the 5 terms, we have only included the four transitions 
indicated in Fig. [2] (i.e., we have neglected the infiuence of the infrared transitions 3d^D- 
3p^P and 3p'^P-3s'^S). The Einstein coefficients for the four included transitions shown in 
Table [3] and the en ergy of the levels s hown in Fig. [2] have b een obtained from the NIST 



databas^ (see also Wiese et al. 1966 : Drake fc Martin 1998). Table [3] also indicates the 



value of the critical magnetic field strength for the operation of the upper- level Hanle effect, 
which results from equating the Zeeman splitting of the level with its natural width: 



^critical ~ 1.137 X 10-V(tlifef7L), (16) 



■^http : //physics .nist .gov/PhysRef Data/ASD/ index . html 
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where tufe is the level's lifetime in seconds, is its Lande factor and -Bcriticai is given in 
gauss. Note that for obtaining the B^^^^^J^^ values of Table [3] we have used tiife~l/v4„/, while 
an estimation of the critical magnetic field strength for the operation of the lower-level Hanle 
effect requires using ~ ^/{Bi^Jq) (which for the metastable lower-level of the He I 10830 



A multiplet gives -B, 



lower 
critical 



0.1 G). 



We point out that in our modeling we are not explicitly taking into account the radiative 
mechanism that is thought to be responsible of the overpopulation of the triplet levels of He I 
required to produce the absorption or emission features observed in the spectral lines of such 
two multiplets -that is, ionizations from the singlet states of He I caused by EUV ionizing 
radiation coming downwards from the corona fo llowed by recombinations towards both the 
singlet and triplet states (e.g.. lAvrett et al.lll994l ). The fact that most of the ionizations take 



place from the singlet states suggests that the atomic polarization of the triplet states should 
be little affected by such EUV coronal irradiation. This expectation is reinforced by the fact 
that the number of photoionizations per unit volume and time from the triplet levels of He I i s 
way smaller than the number of bound-bound transitions (see Fig. 7 in lCenteno et al.ll2008l ). 
which suggests that the atomic polarization of the J-levels of the 10830 A and D3 multiplets 
is indeed dominated by optical pumping in the line transitions themselves. Following our 
approach, the key mechanism responsible of the absorption or emission observed in such 
lines of neutral helium, be it the above-mentioned photoionizatio n- recombination raechan ism 
and/or collisional excitation in regions with T> 20000 K (cf., lAndretta fc Joned 119971 ). is 
unimportant in our forward modeling and inversion approach. The reason is that they are 
implicitly accounted for via the definition of the optical thickness of the slab which, being a 
free parameter i n our modeling, is use d to fit the observed intensity profiles. This point has 
been clarified by lCenteno et al.l (120081 ). 



2.6. The incident radiation field 

As mentioned before, we consider a slab of He I atoms anisotropically illuminated from 
below by the photospheric continuum radiation, which we assume to have axial symmetry 
around the solar local vertical direction. Since the illumination conditions are assumed to be 
known a priori, the radiation field tensors can be calculated directly from the given incident 
radiation. 

For symmetry reasons, it is advantageous to calculate the statistical tensors of the 
radiation field in a reference frame in which the Z-axis is along the vertical of the atmosphere. 
Since the incoming radiation is assumed to be unpolarized, the only non-zero spherical tensor 
components of the radiation field are Jq and Jq. They quantify the mean intensity and the 
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"degree of anisotropy" of the radiation field, respectively. For the case of our plane-parallel 
slab model, their expressions reduce to: 

where fi = cos 6 is the cosine of the heliocentric angle 6. 

It is convenient to parameterize the radiation field in terms of the number of photons 
per mode n and the anisotropy factor w: 

where c is the light speed, h is the Planck constant and u is the frequency of the transition. 
The anisotropy factor fulfills — l/2<w<l. It reaches w = 1 when the radiation is unidirec- 
tional along the polar axis of the reference system, while w = —1/2 when the radiation field 
is azimuth-independent and confined to the plane perpendicular to the quantization axis. 
We have calculated the tensors of the radiation field by using the center-to-limb variation 



and th e wavelength dependence of the solar continuum radiation field tabulated by [Pierce 



(I2OOOI ). Since the He I lines originate in the outer regions of the solar /stellar atmosphere, it 
is necessary to take into account the geometrical effect produced by the non-negligible height 
h above the surface of the star. To this en d, we have applied the strategy outlined in §12.3 



of iLandi Degl'Innocenti fc Landolfil (120041 ). Figure H] shows the sensitivity of w and n to the 
height in the solar atmosphere at which our slab is assumed to be located. The number of 
photons per mode decreases, while the anisotropy factor rapidly increases. This effect is of 
a purely geometrical nature. 

It is imp ortant to note that the theory of spect ral line polarization presented in the 



monograph of iLandi Degl'Innocenti &: Landolfil (120041 ) is only valid if the radiation field il- 



luminating the atomic system is spectrally fiat (independent of frequency) over a frequency 
interval larger than the Bohr frequencies connecting the levels that present quantum coher- 
ences. Fortunately, this is the case for all the lines included in the atomic model of Fig. [2l 
except for the 3p'^P-2s^S transition at 3888.6 A which is situated in a quite crowded region 
of the spectrum. Fortunately, this line is of secondary importance in setting the statistical 



equilibrium, and according to lLandi Degl'Innocenti fc Landolfil (120041 ) the statistical tensors 



can be calculated by simply reducing the photospheric continuum radiation intensity at that 
wavelength by a factor 5. In any case, its inclusion has a rather negligible impact on the 
final results. 



3. The forward modeling code 



By forward modeling we mean the calculation of the emergent Stokes profiles for given 
values of the height h at which the slab is located above the visible solar "surface", of the 
slab's optical thickness and of the magnetic field vector. We have written this option of 
our computer program in a way such that it performs the calculation at various levels of 
realism. We point out that similar type of forward-modeling calculations can be carried out 
also with some of the computer programs mentioned in §1, which are likewise based on the 
density-matrix theory of spectral line polarization. The main difference with ours is that we 
have taken into account radiative transfer effects (without neglecting the magneto-optical 
terms of the Stokes-vector transfer equation) and that we have developed a user-friendly 
interface to facilitate the performing of numerical experiments (see Fig. [5]). In what follows 
we first describe the various options of our computer program and then show some illustrative 
examples of possible applications. 



3.1. Description of the computing options 

The solution of the statistical equilibrium equations, using the radiation field tensors 
calculated from the given illumination conditions, provides the mutipolar components of the 
atomic density matrix. The computer program calculates such Pq ( J, J ) elements in both the 
magnetic field reference frame (if a deterministic magnetic field is present) and in a reference 
system where the quantization axis lies along the solar local vertical direction (hereafter, 
"vertical frame"). This option of the forward modeling code helps to understand what's 
going on at the atomic level when an atomic system is subjected to anisotropic radiative 
pumping processes in the absence and in the presence of a magnetic field. Some illustrative 
examples for the J-levels involved in the 10830 A and D3 transitions are shown in §3.2 below. 

The values of the multipolar components of the density matrix is all we need for cal- 
culating the coefficients of the emission vector and of the propagation matrix, which enter 
the Stokes vector transfer equation ([1]). Our forward modeling code can solve this equation 
at the various levels of approximation explained in §2.11 Moreover, it can also compute the 
emergent Stokes profiles calculating the wavelength positions of the vr and a components 
assuming the linear Zeeman-effect regime (instead of using the general Paschen-Back effect 
theory), incorporating or discarding the influence of atomic level polarization. 
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3.2. Atomic level polarization in two reference systems 

Figure [6] shows an illustrative example of the population imbalances [Po{J) and pl{J)] 
and quantum coherences [pq{J) and Pq{J), with Q^O] induced by optical pumping pro- 
cesses in the levels of the He I 10830 A multiplet that can carry atomic polarization (i.e., 
the lower level, with J; = 1, and the upper levels with J„ = 2 and J„ = 1). Like in 



Trujillo Bueno fc Asensio Ramod (120071 ). we assume a slab with Arred = 1, where the label 
"red" indicates that the optical thickness of the slab along its normal direction is measured 
at the frequency where the peak of the red blended component of the 10830 A multiplet is 
located. The slab is assumed to be at a height of only 3 arcseconds and in the presence of a 
horizontal magnetic field whose strength we can vary at will. 

Consider first the results for Ji = 1 and i? = G in a reference system with the 
quantization axis along the local vertical. Since the pumping radiation is assumed to be 
unpolarized and with axial symmetry around the vertical, for 5 ^ G we see only population 
imbalances of the form Pq{J)- Note that in the absence of magnetic fields the Pq{Ji) atomic 
orientation value is zero because there is no net circular polarization in the incident radiation 
field. The bottom right panel of Fig. [6]shows that, for S = G, we have lower-level quantum 
coherences of the form Pq{J) (with Q^O) in the magnetic reference system, whose Z-axis 
is inclined with respect to the symmetry axis of the pumping radiation field. In fact, for 
given illumination conditions in the absence of a magnetic field, whether or not we have 
such quantum coherences depends only on the reference system. In each of the two bottom 
panels it is shown what happens with the population imbalances and the coherences of the 
Ji = 1 lower level as the strength of the assumed horizontal field is increased. Consider, for 
instance, the right panel results corresponding to the magnetic field reference frame. Note 
that the lower-level Hanle effect starts to operate for field strengths as low as 0.01 G, and 
that for field intensities of the order or larger than 1 G all the lower-level coherences have 
been relaxed. This happens because such a lower level is metastable, which implies that 
its critical Hanle field is only 0.1 G (see Eq. [T6|) . It is also of interest to point out that if 
the J-levels of our multiterm atomic model were isolated levels then the pl{J) population 
imbalances would be constant in the magnetic field reference frame, and the pl{J) atomic 
orientation value would be zero. However, since the J-levels suffer crossings and repulsions 
non-zero p^jJ) values appe ar through the alignment-to-orientation conversion mechanism (cf. 



Landi Degrinnocentilll982l ). which for the 10830 A triplet has however a negligible impact on 
the emergent Stokes V profiles. In addition, as shown in the bottom right panel, the pl{Ji) 
alignment coefficient itself starts to be modified as soon as the field strength is sensibly larger 
than 100 G. Although this lower level does not suffer any crossings with the other J-levels of 
the model atom of Fig. [21 its atomic polarization is modified because it sensitively depends 
on that of the upper levels of the 10830 A multiplet. 
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Consider now the results for the upper levels with = 1 (central panels of Fig. [6]) 
and Ju = 2 (top panels). Obviously, we only see population imbalances of the form pl{J) 
at zero field in the vertical frame. Around B = 10~^ G, the density matrix elements start 
to be affected by the magnetic field. This modification is due to the feedback effect that 
the alteration of the lower-level polarization has on the upper levels. It is also possible to 
note in the central and top panels of Fig. [6] the action of the upper level Hanle effect on the 
Ju = 1 and Ju = 2 levels. In fact, there is a hint of a small plateau just above 0.1 G, as 
expected from the fact that the critical upper-level Hanle field value is 0.8 G for the 10830 
A multiplet (see Table [3]). Similar features can be seen in the corresponding right panels of 
Fig. [6l Probably, the most notable conclusion to highlight from this figure is that between 
approximately 10 and 100 G we only have population imbalances of the form Pq{J) in the 
magnetic field reference frame (i.e., we can consider that between 10 and 100 G the He I 
10830 A multiplet is in the saturation regime of the Hanle effect, where the coherences are 
negligible and the atomic alignment values of the lower and upper levels are insensitive to 
the strength of the magnetic field). 

Finally, in Fig. [7]we show similar results but for the upper levels of the He I D3 multiplet. 
The situation now is much more complicated, as evidenced by the fact that it is impossible 
to find a range of solar magnetic field strength values between which the coherences are zero 
in the magnetic field reference frame and, at the same time, the population imbalances are 
insensitive to the magnetic strength. 



3.3. How to investigate empirically the possibility of magnetic canopies in the 

quiet solar chromosphere? 

In the presence of an inclined magnetic field forward scatteri ng processes can pro duce 
linear polarization signatures in spectral lines (e.g., the review by lTrujillo Buendl200ll ). In 
this geometry, the polarization signal is created by the Hanle effect, a physical phenomenon 
that has been clearly de monstrated via spectropolar imetry of solar coronal filaments in the 
He I 10830 A multiplet JTruiillo Bueno et al.l[2002b[ ). 



Fig. [8] shows theoretical examples of the emergent fractional linear polarization in the 
lines of the He I 10830 A multiplet assuming a constant-property slab of helium atoms 
permeated by a horizontal magnetic field of 10 G. As expected, the smaller the optical thick- 
ness of the assumed plasma structure the smalle r the fractional polarization amplitude. I n 



princ iple, the Tenerife Infrared Polarimeter (TIP; [Martinez Pillet et al.lll999l : ICoUados et al. 



20071 ) mounted on the Vacuum Tower Telescope (VTT) of the Observatorio del Teide allows 
the detection of very low amplitude polarization signals, such as those corresponding to the 
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Arrcd = 0.1 case of Fig. [HI However, in order to be able to achieve this goal without having 
to sacrifice the spatial and/or temporal resolution we need a larger aperture solar telescope. 

We consider now the question of whether it is safer to interpret disk-center observations 
or to opt for a different scattering geometry. To this end, we have investigated how is the 
variation of Q/I and U/I at the central wavelengths of the red and blue components of the 
10830 A multiplet for different inclinations {6b) of the magnetic field vector and for different 
line-of-sights, assuming Axrcd = 0.1 and 5 = 10 G (which implies that we are very close 
to the saturation regime of the upper level Hanle effect). Fig. [9] shows only the case of 
magnetic field vectors with azimuth xb = 90° (i.e., contained in the Z-Y plane of Fig. [1]) 
and for line-of-sights with x = (i-^-, contained in the Z-X plane of Fig. [1]). It is interesting 
to note that for the case of a horizontal magnetic field (i.e., Ob = 90°, which implies that the 
magnetic field vector forms always an angle of 90° with respect to any of the line-of-sights 
of Fig. [9]) Stokes U = and the Stokes Q amplitudes of the emergent spectral fine radiation 
are identical for all such line-of-sights. This is easy to understand by using the following 
approximate expressions for eg/e/ and ?7q/?7/, which for the case of a deterministic magnetic 
field provide a suitable approxima tion if we are in the saturation regime of the upper-level 



Hanle effect ( iTrujillo Buendl2003l ): 

(/x|-l)WK2(J„)]^, (20) 
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(/x^-1)ZK^(J,)]b, (21) 

where [ctqIb = [PqIb/Pq quantifies the fractional atomic ahgnment in the magnetic field 
reference frame, while W and Z are numerical coefficients that depend on the Ji and Ju 
valuefl It is very important to note that in Eqs. (!20|) and (12T1) /is is the cosine of the angle 
between the magnetic field vector and the line of sight. Note also that in these expressions 
for eq/ej and riq/rij, that are valid in the magnetic field reference frame, we have chosen 
the positive reference direction for Stokes Q parallel to the projection of the magnetic field 
ve ctor onto the plane perpendicular to th e LOS, while in the similar Eqs. (16) and (17) 



of ITrujillo Bueno &: Asensio RamosI (120071 ) we chose it along the perpendicular direction. 
In both cases, we have eu = t]u = (if we are really in the above-mentioned Hanle-effect 
saturation regime). Obviously, the observed Q/I amplitude depends on the fractional atomic 



(2) (2) (2) 

Actually, W = w j and Z ~ w j^ j , with w j y given by Eq. (10.12) of lLandi Degl'Innocenti fc Landolfi 



([20041). For instance, W = and Z = 1 for a line transition with Ji — 1 and Ju — 0, such as that of the blue 



component of the 10830 A triplet. 



alignment of the upper and lower levels of the 10830 A line transitions calculated in the 
magnetic field reference frame, but their values are independent of the LOS. Actually, in the 
weak anisotropy limit they are given by 




(22) 



where [(Jq]v is the fractional atomic alignment in the vertical frame for the zero field case. 
Therefore, the dependency of the emergent Stokes Q signal on the scattering angle is only 
through the factor (/i^^ — 1), with fis = for all the line-of-sights of Fig. [9]if 6*^ = 90°. 

As seen in Fig. [HI for non-horizontal fields there are notable differences between the 
curves corresponding to each LOS, mainly for the cases with an inclination {6b) of the 
magnetic field smaller than the Van-Vleck angle {Oyv = 54.73°, which corresponds to 
cos^lOyv) = 1/3). For the forward-scattering case of a disk-center observation (the fi = 1 
curves of Fig. ^ Stokes ~ and Stokes Q admits only one solution for 6*^ > Oyv (with 
the exception of the well-known 180° ambiguity of the Hanle effect )[§. However, for 6b < 6vv 
we may have two different magnetic field inclinations producing the same Stokes Q values, 
each of them having its corresponding 180° ambiguity of the Hanle effect. As seen in the 
figure, the range of values where two such solutions for Stokes Q are possible depends 
on the /i-value of the LOS. Stokes U is clearly non-zero for /i < 1, but the fact that t/ ~ 
for the case of magnetic field with fixed inclination but with random-azimuth within the 
spatio-temporal resolution element of the observation, leads us to conclude that the best one 
can do for a reliable empirical investigation of the possibility of canopy-like fields in the quiet 
solar chromosphere is to interpret spectropolarimetric observations at the solar disk center, 
such as those considered in §5.41 



Our inversion strategy is based on the minimization of a merit function that quantifies 
how well the Stokes profiles calculated in our atmospheric model reproduce the observed 
Stokes profiles. To this end, we have chosen the standard x^^function, defined as: 



4. 



The Inversion Code 




1 



4 



[^r(A,)-^f-(A,)]^ 



(23) 



^The reason why Stokes U in forward scattering geometry is not exactly zero for the red component is 
because for B=10 G some of the coherences are not completely insignificant. 
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where N\ is the number of wavelength points and <yf{\j) is the variance associated to the 
j-th wavelength point of the i-th Stokes profiles. The minimization algorithm tries to find 
the value of the parameters of our model that lead to synthetic Stokes profiles S^^^"^ with 
the best possible fit to the observations. For our slab model, the number of parameters 
(number of dimensions of the hypersurface) lies between 5 and 7, the maximum value 
corresponding to the optically thick case (see Tabled]). The magnetic field vector (5, 6b 
and Xb)i the thermal velocity (fth) and the macroscopic velocity (fmac) are always required. 
This set of parameters is enough for the case of an optically thin slab. In order to account 
for radiative transfer effects, we need to define the optical depth of the slab along its normal 
direction and at a suitable reference wavelength (e.g., the central wavelength of the red 
blended component for the He I 10830 A multiplet). In addition, we may additionally need 
to include the damping parameter (a) of the Voigt profile if the wings of the observed Stokes 
profiles cannot be fitted using Gaussian line profiles. 

A critical problem in any inversion code is to identify possible degeneracies among 
different parameters of the model. When two or more parameters produce similar effects 
on the emergent Stokes profiles, the inversion algorithm is unable to distinguish between 
them. As a result, the emergent Stokes profiles corresponding to different combinations of 
the model parameters are indistinguishable within the noise level. Concerning this critical 
point, we investigate in detail in §6.21 the possibility of using the observed Stokes profiles 
of the 10830 A triplet to obtain the value of the height h at which the observed plasma 
structure is located. 

Ambiguities in the determination of the model's parameters can also result from the 
presence of degeneracies. However, this type of ambiguities occur only for a finite number 
of combinations of some parameters. Although the problem is complicated, it is possible to 
develop techniques that can help selecting one of the combinations as the most plausible. The 
well-known 180° ambiguity of the Hanle effect adds to the unfamiliar Van V leck ambiguity 
flHousel I1977I : ICasini fc Judg3 llQQol : ICasini et all boosl : iMerenda et aP booel ). The role of 
these ambiguities in the inferred model's parameters will be explored in §6.11 although some 
hints have been already given in §3.31 



It is also important to point out the interest of developing methods capable of providing 
reliable confidence intervals for all the inferre d parameters. In a future de velopment, we plan 
to implement Bayesian inference techniques (lAsensio Ramos et al.l 120071 ). 
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4.1. Levenberg-Marquardt 

Probably, the most well-known proce dure for the mini mization of the x^-function is the 



Levenberg-Marquardt (LM) method (e.g. JPress et al.lll986l ). The minimization strategy uses 
the Hessian method when the parameters are close to the minimum of the x^-function (a 
quadratic form approximately describes this function around the minimum) and the steepest 
descent method when the parameters are far from the minimum. The transition between 
both methods is done in an adaptive manner. Its main drawback (which applies also to 
the majority of the standard numerical methods of function minimization) is that it can 
easily get trapped in local minima of the x^-function. Some alternatives are available to 
overcome this difficulty. The most straightforward but time consuming one is to restart the 
minimization process at different values of the initial parameters. If the obtained minimum 
is systematically the same, the probability that this is the global optimum is high. How- 
ever, when the parameter hypersurface is complicated, this technique does not give any 
confidence on the validity of the global minimum. Other possibilities rely on the application 
of some kind of "inertia" to the method, so that the LM method can overcome such a local 
minimum problem when moving on the parameter hypersurface. Again, these methods do 
not guarantee the success of getting the global minimum of the function. On the contrary, 
the LM method turns out to be one of the fastest and simplest options when the initial 
estimate of the parameters is close to the absolute minimum. 



4.2. Global Optimization techniques 



In order to avoid the possibility of getting trapped in a local minimum of the x^ hy- 
persurface, global optimization methods have to be used. Several opti mization methods 
have been developed to obtain the global minimum of a function (e.g., iHorst fc Pardalos 
19951 ). The majority of them are based on stochastic optimization techniques. The essential 
philosophy of these methods is to sample efficiently the whole space of parameters to find 
the global minimum of a given function. One of the most promising methods is genetic 
optimization (inspired by the fact of biological evolution), in which the parameters of the 
merit function are encoded in a gene. Although no mathematical proof of the convergence 
properties of these alg orithms ex i sts, r ecent advances suggest that the probability of con- 
vergence is very high (IGutowskyl 12004 ) . Actually, they perform quite well in practice for 
the optimization of very hard problem s. In solar s pectr opolarimetry, genetic optimization 
methods have been recently applied by iLagg et al.l (120041 ) to the inversion of Stokes profiles 
induced by the Zeeman effect in the He I 10830 A triplet, neglecting the infiuence of atomic 
level polarization. The main disadvantage is that the computing time needed to reach con- 
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vergence increases dramatically (by a factor ~10 with respect to standard methods based 
on the gradient descent like LM). 



^Another different approach is based on deterministic algorithms (IHorst &: Pardalos 

19951 ). Typically, these algorithms rely on a strong mathematical basis, so th at their con- 



verge nce properties are well known. We have chosen the DIRECT algorithm (jjones et al. 



19931 ). whose name derives from one of its main features: dividing rectangles. The idea is to 



recursively sample parts of the space of parameters, improving in each iteration the location 
of the part of the space where the global minimum is potentially located. The d ecision algo 



rithm is based on the assumption that the function is Lipschitz continuous (see iJones et al. 



19931 . for details). The method works very well in practice and can indeed find the minimum 



in functions that do not fulfill the condition of Lipschitz continuity. The reason is that the 
DIRECT algorithm does not require the explicit calculation of the Lipschitz constant but it 
uses all possible values of such a constant to determine if a region of the parameter sp ace 
should be broken into subregions because of its potential interest (see lJones et al.lll993l . for 
details). A schematic illustration of the subdivision process for a function of two parameters 
is shown in Fig. [101 



4.3. Convergence 



Taking into account that the dimension of our space of parameters is between 5 and 7, 
it seems unreasonable to use an algorithm like DIRECT to obtain a precise determination of 
the values of the model's parameters at the global minimum. The reason is that the precision 
in the values of the parameters decreases with the size of the hyperrectangles. Therefore, we 
would need to perform many divisions to end up with a reasonable precision. What we do is 
to let the DIRECT algorithm locate the global minimum in a region whose hypervolume is 
V. This hypervolume is obtained as the product of the length dt of each dimension associated 
with each of the parameters: 



TV 



V = l[di 



(24) 



When the hypervolume decreases by a factor / after the DIRECT algorithm has discarded 
some of the hyperrectangles, its size along each dimension is approximately decreased by a 
factor f^^^. In order to end up with a small region where the global minimum is located, 
many subdivisions are necessary, thus requiring many function evaluations. For this reason, 
it has been observed that although the DIRECT algorithm rapidly finds the region where 
the global mini r num is located, its local conver g ence properties are ra ther poor (see, e.g.. 
Cox et al.l I2OOII : iBartholomew-Biggs et al.l l2002l : iLjungberg et al.l l2004j . for applications in 



- 25 - 



the extremely hard problems of the design of high-speed civil transport, aircrafts and bio in- 
formatics). In summary, the DIRECT method is an ideal candidate for its application as an 
estimator of the region where the global minimum is located, but not for determining it. 

The most time consuming part of any optimization procedure is the evaluation of the 
merit function. The DIRECT algorithm needs only a reduced number of evaluations of the 
merit function to find the region where the global minimum is located. For this reason, we 
have chosen it as the initialization part of the LM method. Since the initialization point 
is close to the global minimum, the LM method, thanks to its quadratic behavior, rapidly 
converges to the minimum. 



4.4. Stopping criterium 

A critical and fundamental problem in the optimization of functions (either local or 
global) is to identify when the method has converged to the solution. We have used two 
stopping criteria for the DIRECT algorithm. The first one is stopping when the ratio between 
the hypervolume where the global minimum is located and the original hypervolume is 
smaller than a given threshold. This method has been chosen when using the DIRECT 
algorithm as an initi alization for the L M method, giving very good results. The other good 
option, suggested by iJones et al. is to stop after a fixed number of evaluations of the 

merit function. 

Since the intensity profile is not very sensitive to the presence of a magnetic field (at 
least for magnetic field strengths of the order of or smaller than 1000 G), we have decided to 
estimate the optical thickness of the slab, the thermal and the macroscopic velocity of the 
plasma and the damping constant by using only the Stokes / profile, and then to determine 
the magnetic field vector by using the polarization profiles. The full inversion scheme, shown 
schematically in Table [2l begins by applying the DIRECT method to obtain a first estimation 
of the indicated four parameters by using only Stokes /. Afterwards, some LM iterations are 
carried out to refine the initial values of the model's parameters obtained in the previous step. 
Once the LM method has converged, the inferred values of fth, v^ac (together with a and 
At, when these are parameters of the model) are kept fixed in the next steps, in which the 
DIRECT method is used again for obtaining an initial approximation of the magnetic field 
vector {B,6b,Xb)- According to our experience, the first estimate of the magnetic field vector 
given by the DIRECT algorithm is typically very close to the final solution. Nevertheless, 
some iterations of the LM method are performed to refine the value of the magnetic field 
strength, inclination and azimuth. In any case, although we have found very good results 
with this procedure, the specific inversion scheme is fully configurable and can be tuned for 
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specific problems. 



5. Applications 

The main aim of this section is to illustrate the application of our inversion code to 
some selected spectropolarimetric observations in the He I 10830 A multiplet, showing that 
it gives results that are in agreement with the published ones. In addition, in §5.41 we show 
a new application aimed at determining the strength and inclination of the magnetic field 
vector in the chromosphere above an internetwork region observed at solar disk center. Note 
that in the following applications Arjed will continue denoting the optical thickness of the 
constant-property slab, along its normal direction, measured at the center of the red blended 
component of the He I 10830 A multiplet. 



5.1. Prominences 



The first application is for the case of solar prominences, which are relatively cool and 
dense plasma structures embedded in the T ~ 10^ K solar corona. In these objects the 
observed Stokes Q and U profiles of the He I 10830 A multiplet are dominated by the 
presence of atomi c level polarizat i on, wh ile the Stokes V profile is dominated by the Zeeman 
effect. Recently, iMerenda et al.l (l2006l ) have shown how to infer the magnetic field that 
confines the plasma of solar prominences via the inversion of the Stokes profiles induced 
by scattering processes and the Hanle and Zeeman effects in the lines of the He I 10830 
A multiplet. They analyzed in detail spectropolarimetric observations of the He I 10830 A 
multiplet in a polar crown prominence and concluded that if the observed prominence was 
located in the plane of the sky the magnetic field had to be relatively strong {B ^ 30 G) 
and inclined by only 25° with respect to the local vertical. 

We have applied our inve rsion code to t h e spe ctropolarimetric observations shown in 
Fig- [m taken from Fig. 9 of iMerenda et al.l (120061 ). We have assumed that the observed 
plasma structure was optically thin and that the prominence is located in the plane of the sky. 
The inversion code was used to infer the value of the thermal velocity Vth, the macroscopic 
velocity shift f mac (to allow for a shift in the wavelength calibration) and the magnetic field 
vector (B,9r,yr)- T he atmospheric height was fixed to ^ = 20", the same value used by 
Merenda et al.l (120061 ) . After the four steps summarized in Table |2l we end up with a thermal 
velocity Vth = 7.97 km s~^, a bulk velocity that is compatible with zero and a magnetic field 
vector characterized hj B = 26.8 G, 6b = 25.5° and xb = 161.0°. These values are in very 
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good agreement with the results of iMerenda et aLl (120061 ). namely B = 26 G, 6b = 25° and 



Xb = 160.5°. Note that since the prominence plasma was assumed to lie in the plane of 
the sky, the following solutions are also valid: 6^15 = 180° — 9b and x*b = ^Xb (i-e., the 
well-known 180° ambiguity of the Hanle effect). 

We point out that the total number of evaluations of the merit function was 132. For 
the inversion of the Stokes profiles corresponding to other points of the field of view, one can 
initialize the inversion using the model's parameters corresponding to the previous point. 
Using a few LM iterations, one should be able to reach the global minimum. In case this 
procedure does not work properly, one should return to the four-steps inversion scheme 
already presented in Table [2l 



5.2. Spicules 



Another interesting problem is the deter r ninati on of the magnetic field vector in solar 
chromospheric spicules. iTrujillo Bueno et al.l (120051 ) interpreted spectropolarimetric obser- 
vations of spicules in the He I 10830 A multiplet and concluded that the magnetic field of 
spicules in quiet regions of the solar chromosphere has a strength of the order of 10 G and 
is inclined by about 35° with respect to the local vertical. Their conclusion that the typical 
magnetic field strength is ~ 10 G required to obtain the longitudinal component of the mag- 
netic field ve ctor via some careful m easurements of the Stokes V profiles, such as that shown 
in Fig. 13 of iTrujillo Buend (120051 ). This was needed because for field strengths larger than 
a few gauss the He I 10830 A multiplet enters the saturation regime of the upper-level Hanle 
effect and the observed Stokes Q and U profiles provide only information on the orientation 
of the magnetic field vector. In fact, the application of our inversion code to the observed 
Stokes profiles shown in Fig. [12] (where the Stokes V profile is at the noise level), assum- 
ing that the spicular material is located in the plane of the sky, provides several different 
magnetic field vectors that lead to equally reliable fit s. One of them, given by B = 10 G, 
6b = 37° and Xb = 172° is similar to the one chosen by lTrujillo Bueno et al.l (120051 ). Another 
possible fit is the one illustrated in Fig. [121 which corresponds to B = 2.6 G, 6b = 37° and 
XB = 35°. 



Fig. [13] gives the values of the function for all possible magnetic field inclinations 
and azimuths corresponding to the cases i? = 10 G (left panel) and B = 2.6 G (right panel). 
In each of these panels we have indicated with white dots and numbers the four solutions 
that correspond to equally good best fits to the observed Stokes profiles of Fig. [12] The 
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pair of solutions 1 and 2 correspond to the Van-Vleck ambiguitjo- The same apphes to 
the 1' and 2' solutions. On the contrary, the pair of solutions 1 and 1' or the 2 and 2' are 
not strictly equivalent. The inversion code considers such pairs of solutions as equivalent 
because the observed Stokes V profile is at the noise level and it is not able to differentiate 
between the two cases. Note that these pairs of solutions give exactly the same Stokes Q 
and U profiles, but their corresponding Stokes V profiles have opposite signs. Concerning 
each pair of solutions in Fig. [131 if is possible to verify that the projections on the plane of 
the sky of the magnetic fields corresponding to solutions 1 and 2 form an angle close to to 
90°, which is typical of the Van-Vleck ambiguity. The same happens for the magnetic fields 
corresponding to solutions 1' and 2'. This holds for both cases, i? = 10 G and B = 2.6 G. 
As pointed out above, when the observed Stokes V signal is very small, it is very hard (or 
impossible) to differentiate between the two possibihties having azimuths xb and 180° — Xb- 
The magnetic field vectors 1 and 1' (or those corresponding to the 2 and 2' solutions of 
Fig. [T3|l have the same projection on the line of sight, except for a sign change. Therefore, 
the detection of Stokes V turns out to be furidame ntal to determine which is the correct 
one (jTrujillo Bueno et al.ll2005l : iMerenda et al.ll2006l ). Apart from the considered solutions, 
which are restricted to the interval 0° < < 180°, one has also to take into account the 
well-known ambiguity of the Hanle effect, which applies when the emitting plasma is located 
in the plane of the sky. In this case, we have to add to the possible set of solutions all the 
combinations fulfilling 0^ = 180° — 6b and x*b = ~Xb since both pairs produce the same 
Stokes profiles|l|. 



5.3. Filaments 



We have also considered the inversion of the Stokes profiles presented in lTrujillo Bueno et al 



(l2002bl ). which were observed in a solar coronal filament at the solar disk center. Such pro- 
files, which are reproduced in Fig. [TH were used by those authors to demonstrate the 
presence of atomic polarization in the lower level of the 10830 A multiplet and that the 



^Infor mation on this ambigui ty, typical of the Hanle-efFect saturation regime, can be found in lCasini et al 



(|2005l) . in lMerenda et all (|2006l ) and in fgH] below 



^We point out that this ambiguity of the Hanle effect applies only to some particular scattering geometries 
(i.e., to those of 90° and of forward scattering). If the observed plasma structure is not located in the plane 
of the sky (which implies scattering processes at an angle different from 90°), or if it is not located at the 
solar disk center, then one has a sort of quasi-degeneracy which can disap pear when the angle 9 of Fig. 1 is 
consid erably different from 90° or from 0°. This fact has been exploited bv lLandi Degl'Innocenti fc Bommier 
(|l993l ) to propose a method for removing the azimuth ambiguity intrinsically present in vector magnetograms. 
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Hanle effect due to an inclined magnetic field creates linearly polarized radiation in forward 
scattering geometry. Note that the Stokes Q, U and V profiles are normalized to the maxi- 
mum depression in Stokes / (which is 0.4 Ic, appr oximately) . The appli c ation o f our inversion 
code using h = 40" confirms the conclusions of iTrujillo Bueno et al.l (l2002b( ). yielding the 
following values for the model's parameters: At = 0.86, Vth = 6.6 km s~^, a = 0.19 and a 
magnetic field vector characterized by S = 18 G and 6b = 105°. 



5.4. Inter-network regions 

As discussed in §3.31 the investigation of the possibility of horizontal magnetic canopies 
in the quiet solar chromosphere above internetwork regions is feasible with TIP, especially 
when interpreting measurements of the polarization of the He I 10830 A multiplet in forward 
scattering at the solar disk centeij^. We present in Fig. [15] an observation carried out with 
TIP very close to the solar disk center (/x = 0.98). The slit was crossing an enhanced network 
region of circular shape. A time series of 50 steps with an integration time of 3 seconds was 
performed. The resulting polarimetric sensitivity after averaging over the 50 time steps and 
along a small spatial interval within the observed internetwork region is close to 6xl0~^ in 
units of the continuum intensity. According to the results of the right panel of Fig. [HI this is 
sufficient for detecting the linear polarization signal of a horizontal magnetic field provided 
the optical depth at the wavelength of the red component of the 10830 A multiplet is larger 
than ~ 0.01. The reference system has been rotated until Stokes U is minimized. Since the 
inferred magnetic field strength implies that the He I 10830 A is in the saturation regime of 
the upper-level Hanle effect, the resulting reference direction for Stokes Q lies either along the 
projection of the magnetic field vector on the solar disk, or along the direction perpendicular 
to such a projection. We have applied our inversion code to the above-mentioned observed 
profiles assuming a slab located at a height of 3 arcsec and we have obtained the following 
results: At = 0.19, Vtt = 9.2 km s~^, a = 0.62 and a magnetic field vector characterized by 
B = 35 G, 9b = 21° and xb = 0°. However, other possible solutions can be found with a 
similar goodness of the fit (e.g., B = 47 G, 9b = 47° and xb = 0°). These results obtained 
from a solar disk center observation suggest the presence of magnetic fields inclined by no 
more than 50° in the observed quiet Sun chromospheric region. 



For a preliminary interpretation of an observation of a quiet region located at fi ~ 0.5 see 



Lagg 



POOTI) . 
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5.5. Emerging magnetic flux regions 



As pointed out by iTrujillo Bueno &: Asensio Ramoa (120071 ). the modeling of tlie emer- 
gent Stokes Q and U profiles of the He I 10830 A multiplet should be done by taking into 
account the possible presence of atomic level polarization, even for magnetic field strengths 
as large as 1000 G. An example of a spectropolarimetric observation of an emerging magnetic 
flux region is shown by the circles o f Fig . [T6l The solid lines show the best theoretical fit 
to these observations of iLagg et al.l (120041 ). Here, in addition to the Zeeman effect, we took 
into account the influence of atomic level polarization. The dotted lines neglect the atomic 
level polarization that is induced by anisotropic radiation pumping in the solar atmosphere. 
Our results indicate the presence of atomic level polarization in a relatively strong field re- 
gion (~1000 G). However, it may be tranquilizing to point out that both inversions of the 
observed profiles yield, at least for this case, a similar magnetic field vector, in spite of the 
fact that the corresponding theoretical fit is much better for the case that includes atomic 
level polarization. 



6. Ambiguity and degeneracies 

In the previous subsections we have shown how our inversion code can be used for 
recovering the parameters of the assumed slab atmospheric model from the Stokes profiles 
observed in different solar plasma structures. The aim of this section is to discuss the Van 
Vleck ambiguity and to investigate whether we can infer the height of the observed plasma 
structure directly through the inversion approach. 



6.1. Van Vleck Ambiguity 



In general, the solution to any inversion problem is not unique -that is, it is of- 
ten possible to detect seve ral solutions which are compatible with the observations (e.g., 
Asensio Ramos et al.ll2007l ). Some of the unicity problems are associated with the physics of 
the polarization phenomena (e.g., the 180° ambiguity of the Hanle effect for plasma struc- 
tures located in the plane of the sky or the Van Vleck ambiguity). However, as seen in §5.2[ 
in addition to this type of ambiguities, other degeneracies can appear because of the presence 
of noise in the observed Stokes profiles. 

The Van Vleck ambiguity occurs only for some combinations of the inclinations and 
azimuths. Moreov er, it occurs raainly in the saturation regime of the Hanle effect. For 
example. Fig. 6 of iMerenda et al.l (120061 ) shows the region of parameters for which the Van 
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Vleck ambiguity occurs in the Hanle-effect saturation regime of the He I 10830 A triplet. 
Since two different magnetic field vectors give rise to exactly the same emergent Stokes 
profiles, it is impossible to distinguish between them using only the 10830 A multiplet (or 
four solutions, if the 180° ambiguity of the Hanle effect also applies). However, if more 
information is introduced in the inversion procedure (for instance, by using simultaneous 
observations in the 10830 A and D3 multiplets), it might be possible to distinguish between 
the two possible solutions. 

Unfortunately, it is not easy to determine the ran ge of parameters in which we may 



have the Van Vleck ambiguity. One possibility (used by lMerenda et al.ll2006l ) is to consider 
the theoretical Hanle diagram of the red line of the He I 10830 A multiplet and detect if 
the observed profiles fall in the ambiguity region. We propose another method based on 
the DIRECT algorithm implemented in our inversion code. The DIRECT algorithm can 
rapidly detect regions of the space of parameters where the global minimum may be located. 
Therefore, we can take advantage of this property to identify the two (or more) points in 
the space of parameters {9b, Xb) that produce the same emergent Stokes profiles for a given 
magnetic field strength. 

To this end, we have calculated the synthetic emergent Stokes profiles of the He I 10830 
A line from an optically thin prominence, located in the plane of the sky, with Vth = 8 
km s~^, h = 20", B = 25 G, 6b = 40° and xb = 19°. According to the Hanle diagram 



shown by iMerenda et al.l (120061 ) . these profiles are indistinguishable from the ones given by 
the combination B = 22 G, 6b = 100° and xb = 46°. To these two combinations, we have 
to add those corresponding to the 180° ambiguity: {B = 25 G, 6b = 140°, xb = —19°) 
and {B = 22 G, 6b = 80°, xb = —46°). Using the standard four steps inversion procedure 
explained in ^ the global minimum is rapidly located at position B = 22 G, 6b = 100° 
and Xb = 46°. Keeping fixed the value of all the thermodjTiamical properties and the 
field strength, the DIRECT algorithm is used to recover the inclination and azimuth of the 
magnetic field vector. We show in Fig[T7]the position in the {6b, Xb) space of parameters of 
the N evaluations performed by the DIRECT method. It has been possible to detect the two 
combinations of parameters that give the same emergent Stokes profiles, as stated above. 
With only N = 100 evaluations of the merit function, the DIRECT algorithm has located 
and refined the position of the global minimum. It has also identified the second possible 
solution. When the number of function evaluations increases (even with only = 200), the 
DIRECT algorithm rapidly locates the two minima. For > 200, we face a degradation in 
the convergence rate as discussed in §4.31 

As already discussed, an interesting property of the DIRECT method is that no hyper- 
rectangle is ever discarded from the search. Therefore, a rectangle that in one iteration is 
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not considered to be potentially interesting, can be chosen for division in a later iteratioio. 
This behavior is shown in Fig. [T71 When = 100, only a part of the space of parameters 
has been sampled, with clear gaps for inclinations above 110°. In spite of these gaps, the two 
global minima have been already found. However, when the number of function evaluation 
is increased, the numerical scheme finally evaluates the function in those regions with the 
aim of discarding the presence of additional global minima. 



6.2. Can we infer the height of the observed plasma structure? 

In this section we briefly discuss the possibility of determining the height at which the 
He I atoms are located by only using the information contained in the Stokes profiles of the 
10830 A multiplet. For simplicity, we consider first the optically thin limit, the case of off- 
limb observations (i.e., 90° scattering geometry) and a magnetic field with a fixed azimuth 
and strength. The synthetic profiles correspond to the case v^h = 8 km s~^ and h = 20", 
with the magnetic field vector given by 5 = 25 G, 6b = 40° and xb = 19°. The aim of this 
experiment is to infer the inclination 6b and height h from synthetic Stokes profiles without 
noise. The positions where the merit function has been evaluated by the DIRECT algorithm 
are presented in the upper panels and in the bottom left panel of Fig. [181 The surface 
is shown in the bottom right panel of Fig. [TSl The presence of the vertical strip where the 
minimum is located makes it very difficult to converge to the minimum using gradient-based 
methods like the LM method. The derivatives cannot be correctly approximated when the 
function has a large variation in such a small region of the space of parameters. This shallow 
strip is produced by the quasi- degeneracy of the problem in both parameters. An infinity of 
combinations of both parameters give Stokes profiles that can approximately reproduce the 
observations. The difference in the merit function between these spurious cases and the 
exact one is very small. Two reasons produce such a behavior. On the one hand, the linear 
polarization signal is enhanced when the height is increased because the anisotropy of the 
radiation field increases (see Fig. HJ right panel). On the other hand, the Hanle effect turns 
out to be particularly efficient in reducing the atomic polarization when the magnetic field 
is significantly inclined with respect to the symmetry axis of the radiation field (the vertical 
direction). In a realistic case, the problem is much more complicated due to the presence of 
other additional parameters and the noise contamination. 

When the observed structure is off the limb, imaging techniques can be used to esti- 



^T his proves to be fundamental to demonstrate that the global minimum will always be found ( Jones at al. 
19931) . 
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mate h. On the contrary, the case of on-disk observations is much more comphcated since no 
straightforward technique for estimating the height is available. One possibility is to follow 
the observed active region until it approaches the limb. The height can then be estimated 
if we assume that the height of the plasma structure has not changed between both obser- 
vations. An even less precise procedure is to assume a given h value based on the typical 
height of the solar structure type under study. Obviously, the ideal situation would be the 
one where h could be inferred directly from the observed Stokes profiles. In order to inves- 
tigate this possibility, we have performed an experiment in which the DIRECT method is 
used with disk-center {6 = 0°) synthetic Stokes profiles. The emergent profiles have been 
calculated with Wth = 8 km At = 0.8, h = 20", B = 25 G, 63 = 40° and xb = 19°, tak- 
ing into account the effects of radiative transfer in the slab. We keep fixed all the parameters 
except for the inclination of the magnetic field 6b and the height h. The upper panels and 
the bottom left panel of Figure [12] present the points at which the DIRECT algorithm has 
evaluated the merit function, showing that it is possible to infer the height of the observation 
by only using the Stokes profiles. The shape of the surface is shown in the bottom right 
panel of Fig. [191 In comparison with the off-limb case shown in Fig. [TH the minimum is 
located in a much less complicated region of the surface. The quasi-degeneracy present 
in the off-limb case is not present in the on-disk case. This is associated with the fact that 
the blue component gives no signal in the optically thin liniit, wh ile it does if an inclined 



field is present for the disk center case (iTrujillo Bueno et al.ll2002bl ). 



Interestingly, if one wants to infer the magnetic field vector and the height simultane- 
ously from the observations, the code is unable to get a suitable global minimum, even in 
the noiseless case. However, an easily accessible global minimum exists when one of the 
parameters is kept fixed, thus inferring only the following combinations of parameters: {B, 
9b, h), {B, Xb, h) and {Ob, Xb, h). 



7. Conclusions 

The physical interpretation of spectropolarimetric observations of lines of neutral he- 
lium, such as those of the 10830 A and D3 multiplets, represents a very important diagnostic 
window for investigating the dynamical behavior and the magnetic field of plasma struc- 
tures in the solar chromosphere and corona, such as spicules, filaments, regions of emerging 
magnetic fiux, network and internetwork regions, sunspots, fiaring regions, etc. In order 
to facilitate this type of investigations we have developed a powerful forward modeling and 
inversion code that permits either to calculate the emergent spectral line intensity and polar- 
ization for any given magnetic field vector or to infer the dynamical and magnetic properties 
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from the observed Stokes pr ofiles. This diagnostic tool is based on the quantum theory of 



spectral line polarization (see lLandi Degl'Innocenti fc Landolf]||2004j ). which self-consistently 



accounts for the presence of atomic level polarization and the Hanle and Zeeman effects in 
the most general situation of the incomplete Paschen-Back effect regime. It is also of interest 
to mention that the same computer program can be easily applied to other chemical species 
apart from He I (e.g., in order to investigate the magnetic sensitivity of the polarization 
caused by the joint action of the Hanle and Zeeman effects in many other spectral lines of 
diagnostic interest, both in the solar atmosphere and in other astrophysical plasmas). 

The influence of radiative transfer on the emergent spectral line radiation is taken into 
account by solving the Stokes-vector transfer equation in a slab of constant physical proper- 
ties, including the magneto-optical terms of the propagation matrix. Although this "cloud" 
model for the interpretation of polarimetric observations in such lines of He I is suitable 
for inferring the magnetic field vector of plasma structures embedded in the solar chromo- 
sphere and corona, there are several interesting improvements and generalizations on which 
we are presently working on. The first one will be useful for improving the modeling of the 
Stokes profiles observed in low-lying optically-thick plasma structures embedded in the solar 
chromosphere, such as those of active region filaments. It consists in taking into account 
that in optically-thick plasma structures located at low atmospheric heights, the atomic level 
polarization is not going to be necessarily dominated by the anisotropic continuum radia- 
tion coming from the underlying solar photosphere (as we have assumed here), given that 
the radiation field generated by the optically-thick structure itself will tend to reduce the 
anisotropy factor of the true radiation field that pumps the helium at oms of the plasma struc- 



ture under consideration (see lTrujillo Bueno fc Asensio Ramosll2007l ). The second additional 



development consists in considering a Milne-Eddington atmospheric model, but determining 
consistently the height-dependent atomic level polarization induced by the anisotropic ra- 
diation field within the atmosphere model that provides the best fit to the observed Stokes 
profiles. S ince the anisotropy factor is very sensitive to the source-function gradient (e.g.. 
Fig. 4 in iTruiillo Bue^ boOlh the solution of these type of problems in stratified model 



atmospheres r aay be facilitated by the application of efficient iterative schemes, such as 
those used by iManso Sainz fc Trujillo Bueno J2003al lb[) for developing a general multilevel 



radiative transfer program for modeling scattering line polarization and the Hanle effect in 
weakly magnetized stellar atmospheres. 

For the solution of the Stokes inversion problem we have applied an efficient algorithm 
based on global optimization methods, which permits a fast and reliable determination of 
the global minimum and facilitates the determination of the solutions corresponding to the 
unfamiliar Van-Vleck ambiguity. Our inversion approach is based on the application of 
the Levenberg-Marquardt (LM) method for locating the minimum of the merit function 
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that quantifies the goodness of the fit between the observed and synthetic Stokes profiles. 
However, gradient-based methods suffer from convergence problems when the initial value of 
the parameters is not close to the minimum. In order to improve the convergence properties 
of the LM method, we have introduced a novel initialization technique. This method is based 
on the DIRECT algorithm, a deterministic global optimization scheme that performs very 
well. We have shown that a four-steps scheme using the DIRECT method to initialize the 
parameters and the LM method to refine the first estimation close to the minimum leads to 
a very robust technique. 

Our computer program has been developed with the aim of being computationally 
efficient and user-friendly. The relevant equations of the problem result from a general 
and robust theory, so that it is straightforward to treat limiting cases and include or discard 
several physical effects in a very transparent way. It is appropriate for its application to a 
wide variety of problems, from simple Zeeman-dominated Stokes profiles to more complex 
situations in which the influence of atomic level polarization cannot be neglected. The code 
is written in FORTRAN 90, and incorporates a user-friendly front-end based on IDL0 which 
facilitates the execution and analysis of the synthesis and/or inversion calculations (see Fig. 

ED. 

Obviously, our inversion strategy cannot compete in spee d with algorithms bas ed on 



look-up tables, like those applied by ICasini et al.l (120031 ) and iMerenda et al.l (120061 ). At 



present, with a modern portable computer, we need of the order of 1 min. for the inversion 
of the Stokes profiles shown in Fig. [T51 The strength of our approach is that it is very 
general and robust, and very suitable also to investigate the impact of the different physical 
mechanisms and parameters on the retrieved models. Concerning future improvements, we 
think that it would be worthwhile to treat the inversion prob lem within the framework of 



Bayesian inference techniques (see lAsensio Ramos et al.ll2007l . for a first application of such 



techniques to the inference of Milne-Eddington parameters from Stokes profiles induced by 
the Zeeman effect). The aim is to sample the joint posterior probability distribution of the 
parameters of the model once the observation has been taken into account, and to carry out 
marginalizations to infer the probability distribution of each parameter. One of the main 
obstacles to overcome is to determine how to sample efficiently the full posterior probability 
distribution in the complex physical problem that we have investigated in this paper. A 
possible solution could be to rely on machine learning techniques for a fast solution of the 
forward p roblem, someth ing that could be in perfect synergy with Markov Chain Montecarlo 



methods (IMacKavl 120031 ). 
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The reliability of the developments presented in this paper has been demonstrated 
through several model calculations and applications. Of particular interest is the investi- 
gation described in Section 3.4, which aimed at clarifying which is the optimum strategy for 
determining, from He I 10830 A spectropolarimetric observations, whether or not we have 
magnetic canopies with horizontal fields in the quiet solar chromosphere. The results of an 
aplication to an observation of a disk-center internetwork region can be found in §5.4, which 
suggest the presence of magnetic fields inclined by no more than 50° in the observed quiet 
chromospheric region. 

We have also discussed the potential problems that one may encounter. For example, we 
have investigated the presence of degeneracies, paying particular attention to the possibility 
of determining the height of the observed plasma structure from the observed Stokes profiles 
themselves and to demonstrate that the DIRECT method is a very efficient technique for 
detecting the solutions associated to the Van Vleck ambiguity. 

"HAZEL" (an acronym from HAnle and ZEeman Light) is the name we have given to 
our lAC computer program for the synthesis and inversion of Stokes profiles resulting from 
the joint action of the Hanle and Zeeman effects in slabs of finite optical thickness. HAZEL 
will be continuously improved over the years (e.g., with extensions to more complicated 
radiative transfer models), but is now ready for systematic applications to a variety of spec- 
tropolarimetric observations in the spectral lines of the He I 10830 A and D3 multiplets. We 
offer it to the astrophysical community with the hope that it will help researchers to achieve 
new breakthroughs in solar and stellar physics. To get a copy, it suffices with making an 
e-mail request to the authors of this paper. 

Acknowledgments We thank Roberto Casini (HAO) for carefully reviewing of our 
paper. Finantial support by the Spanish Ministry of Education and Science through project 
AYA2007-63881 and by the European Commission through the SOLAIRE network (MTRN- 
CT-2006-035484) is gratefully acknowledged. 
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Fig. 1. — The geometry for the scattering event. The Z-axis is placed along the vertical to 
the solar atmosphere. The magnetic field vector, B. is characterized by its modulus -B, the 
inclination angle 9b and the azimuth xb- The linc-of-sight. indicated by the unit vector fl, 
is characterized by the two angles 9 and x- The reference direction for Stokes Q is defined by 
the vector ei on the plane perpendicular to the line-of-sight. This vector makes an angle 7 
with respect to the plane formed by the vertical and the line-of-sight. In the figures showing 
examples of the emergent Stokes profiles, our choice for the positive reference direction for 
Stokes Q is 7 = 90°, unless otherwise stated. For off-fimb observations, we have 9 — 90°, 
while for observations on the solar disk, we have 9 < 90°. Note also that x is generally taken 
to be 0°. 
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Fig. 2. — Model atom of the triplet system of He I used in this investigation. This work 
focuses on the polarization properties of the 10830 A multiplet between the 2p^P and 2s^S 
terms and on the D.s multiplet betwee n the Sd'^D and 2p^P terms. The energy of each J-level 
is taken from lDrake &: Martini (119981 ) and it is given in cm^^ above the fundamental energy 
level (Is^ ^So). Note that the separation between the J-levels pertaining to each term is not 
drawn to scale. 
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Fig. 3. — Energy splitting due to the presence of a magnetic field for the upper levels of the 
10830 A multiplet (left panels) and for the upper levels of the D3 multiplet at 5876 A (right 
panels). Note that the regimes where level crossings and repulsions occur are different for 
the two terms. They can be better identified in the lower panels. The energy of each level 
is referred to the energy of the level with the smallest value of J at zero magnetic field (i.e., 
J = for the upper levels of the 10830 A multiplet and J = 1 for the upper levels of the D3 
multiplet). The energy separation at -B = G is obtained from the information presented 
in Fig. H 
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Fig. 4. — These two quantities (see Eqs. [I9]) characterize the radiation field that produces 
optical pumping processes in the He I atoms. The number of photons per mode n [left 
panel) is proportional to the mean intensity, Jq, while the anisotropy factor w [right panel) 
is proportional to the Jq tensor of the ra diation f ield. Both quantities have been obtained 
using the limb-darkening data tabulated by lPiercd (120001 ). The figures show also the variation 
of these quantities with the atmospheric height at which the slab of helium atoms is assumed 
to be located. Note that, while the number of photons per mode is almost insensitive to h, 
the geometrical effects produce a significant variation on the anisotropy factor. 
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Fig. 5. — Screen dump of the graphical front-end used for the synthesis. 
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Fig. 6. — Variation of the fractional population imbalances, Pq{J)/pq{J), and of the non-zero 
quantum coherences between magnetic sublevels pertaining to a given J-level, Pq {J) / PoiJ) , 
for different values of the strength of a horizontal magnetic field. We show these quantities 
for the Ju = 1 and ,/„ = 2 levels of the upper term of the 10830 A transition (upper and 
middle panel, respectively) and for the J; = 1 level of the lower term of the 10830 A transition 
(lower panel). The left panels show the results for the "vertical" reference frame in which 
the quantization axis is chosen along the symmetry axis of the radiation field, while the right 
panels show the results for the magnetic field reference frame in which the quantization axis 
is chosen along the (horizontal) magnetic field vector. Note that, at B^iO G, only population 
imbalances are present in the vertical reference frame, while both population imbalances 
and coherences are present in the magnetic reference frame. Note that in the magnetic field 
reference frame the quantum coherences are zero and pl{J) is constant for 10 < B < 100 G. 
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Fig. 7. — Same as Fig. [6l but for the J-levels of the upper term of the He I D3 multiplet. 
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Fig. 8. — Emergent fractional linear polarization at the solar disk center in the lines of 
the He I 10830 A multiplet, assuming that a constant-property slab of helium atoms at a 
height of 3 arcseconds is permeated by a horizontal magnetic field of 10 G. The various Q/I 
profiles of the left panel correspond to the slab's optical thickness indicated in the inset, 
calculated at the wavelength of the red blended component. The right panel indicates that 
the fractional polarization amplitudes of the red and blue components increase exponentially 
with the slab's optical thickness, till line saturation sets in. These calculations have been 
carried out using the exact solution of the radiative transfer equation given by Eq. (j5]). 
The positive reference direction for the definition of the Stokes Q parameter lies along the 
horizontal magnetic field vector. 
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Fig. 9. — Variation of the linear polarization at the line center of the blue (left panel) and 
red (right panel) components of the He I 10830 A multiplet for different inclinations of the 
magnetic field vector with respect to the local vertical and for different observing angles. The 
upper panels show Stokes Q while the lower panels show Stokes U , normalized to the value of 
Stokes / at the fine center of each component. The calculations have been obtained assuming 
that a constant-property slab of helium atoms at a height of 3 arcseconds is permeated by a 
magnetic field of 10 G with an inclination 9b with respect to the solar local vertical direction. 
The slab's optical thickness at the wavelength of the red blended component is Arred = 0.1. 
The positive reference direction for Stokes Q is along the projection of the magnetic field 
vector on the solar surface, which makes an angle of 90° with any of the considered line-of- 
sights. 
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Fig. 10. — This figure illustrates the philosophy of the DIRECT method for searching the 
region where the global minimum is located. In this case, we present an illustrative example 
in two dimensions. After the evaluation of the merit function at some selected points inside 
each region, the DIRECT algorithm decides, using the Lipschitz condition, which rectangles 
should be further subdivided in case they are potentially optimal (optimal means that a 
low value of the merit function has been found or that the rectangles are large and a finer 
sampling has to be carried out). This method rapidly finds the region where the minimum 
is located. 
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Fig. 11. — The solar prominence case. Example of the best theoretic al fit obtained with th e 
inversion code to some of the observed Stokes profiles presented in iMerenda et al.l (120061 ). 
The emergent Stokes profiles have been obtained by assuming an optically thin plasma. The 
inferred magnetic field vector is i? = 26.8 G, 6b = 25.5° and xb = 161° and the inferred 
thermal velo city Vth = 7.97 km s~\ All these values are in good agreement with the results 
obtained by iMerenda et al.l (120061 ). The positive reference direction for Stokes Q is the 
parallel to the solar limb. 
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Fig. 12. — The case of solar chromospheric spicules. Example of th e best theoretical fit 



obtain ed with the inversion code to the observed profiles presented in iTrujillo Bueno et al. 



(120051 ) ■ The emergent Stokes profiles have been obtained by taking into account radiative 
transfer effects in a constant-property slab model. The inferred magnetic field vector is 
B = 2.6 G, 9b = 37° and xb = 35°, although an equally good fit is obtained for other 
combinations like i? = 10 G, = 37° and xb = 172° (see the text for more information). 
Furthermore, the inferred thermal velocity is tith = 13.9 km s~^, the slab's optical thickness 
in the red blended component is Arred = 2.54 and the damping parameter is a = 0.22. The 
positive reference direction for Stokes Q is the parallel to the solar limb. 
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Fig. 13. — Values of t he Y^-function for the possi ble fits of the Stokes profiles of the spicular 
material observed by lTrujillo Bueno et al.l (120051 ). after considering several combinations of 
the inclination and azimuth of the magnetic field vector. We fixed the thermal velocity, 
the optical depth of the slab and the damping constant to the following values: fth = 13.9 
kms~^, ATj-cd = 2.54 and a = 0.22. The magnetic field strength is 10 G in the left panel 
while it is 2.6 G in the right panel. Note the presence of several local minima. The white 
dots indicate the position of the global minimum of the x^-function. The solutions 1 and 
2 correspond to the Van-Vleck ambiguity. The same happens for the solutions 1' and 2'. 
Additionally, the solutions 1 and 1' are equivalent for the code even if they do produce a 
sign change in the Stokes V signal because the observed Stokes V profile is at the level of 
noise. The same happens for the solutions 2 and 2'. 
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Fig. 14. — The case of a solar coronal filament. Example of the Stokes profiles of the He I 
10830 A multiplet observed in a corona l filame nt at solar disk center. Our results confirm 
the conclusions of iTrujillo Bueno et al.l ( j2002bl ). since we obtain At = 0.86, Vth = 6.6 km 
s~^, a = 0.19, B = 18 G and 6b = 105°. The positive direction of Stokes Q is parallel to the 
projection of the magnetic field vector on the solar surface. 
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Fig. 15. — The quiet chromosphere case. Example of the Stokes profiles of the He I 10830 
A multiplet observed in a disk center (/j, — 0.98) internetwork region surrounded by an 
enhanced network region of almost circular shape. The observed Stokes profiles selected for 
the inversion correspond to a temporal and spatial average within the internetwork region. 
The solid line presents a fit to the observations which corresponds to a magnetic field vector 
with B = 35 G, 6b = 21° and xb = 0°. The inferred optical depth in the red blended 
component is Arrcd ~ 0.2, while fth = 9.2 km s^^ and a = 0.62. An equally good fit is 
obtained for other field configurations like B = 47 G, 6b = 47° and xb = 0°. The fact that 
with &6b < 54.74° the Stokes Q signal of the red component is positive (after our rotation of 
the reference system to minimize Stokes U) indicates that the reference direction for Stokes 
Q lies along the direction perpendicular to the projection of the magnetic field vector on the 
solar disk (see the ji — l case of Fig. 9). 
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Fig. 16. — The case of an emerging flux region. Example of the Stokes profiles of the 
He I 10830 A multiplet in an emerging flux region lo cated at u = 0.8. The circles present 
the observed Stokes profiles obtained from Fig. 2 of iLagg et al.l (120041 ) . but after rotating 
the reference system by —56° to have the reference direction for the Stokes Q profile along 
the direction perpendicular to the straight line joining the disk center with the observed 
region. For this reason the inversion has been carried out using x = 0° and 7 = 90° (see 
Fig. 1). The solid line presents the best fit obtained with our inversion code, taking into 
account the effect of atomic polarization on the emergent Stokes profiles, in addition to the 
Zeeman effect treated within the framework of the Paschen-Back effect theory. This full 
solution corresponds to a magnetic field vector with B = 1073 G, 6b = 86° and xb = 
170.7°, while Arrcd = 1.07, Vth = 7.24 km s~^ and a = 0.25. The dotted line is the best 
fit obtained when neglecting atornic po larization. These results confirm the conclusion by 
Trujillo Bueno fc Asensio Ramod (120071 ) that the presence of atomic polarization in the levels 
of the He I 10830 A multiplet produces observable signatures on the emergent Stokes profiles 
even for fields as large as 1000 G. We point ou t that there is a typing error in the caption 
of figure 3 of iTrujillo Bueno &: Asensio RamosI ( 120071 ) . since the theoretical Stokes profiles 
in that figure were obtained for a LOS with n = cosd = 0.8 (i.e., not for the /i = 1 forward 
scattering case) and with x = 7 = 0° (see Fig. 1). 
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Fig. 17. — Value of 9b and xb at which the merit function is evaluated for detecting possible 
ambiguities using the DIRECT method. The number of function evaluations is shown on 
the top of each panel. The two global minima seen in the figures result from the Van Vleck 
ambiguity, showing two different magnetic field vectors that produce the same emergent 
Stokes profiles. This figure illustrates the power of the DIRECT method, which allows to 
recover the two values of the magnetic field vector in this ambiguous case with less than 200 
function evaluations. We have shown only half of the space of parameters in order to avoid 
the 180° ambiguity. 
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Fig. 18. — Value of 6b and h at which the merit function is evaluated when using the 
DIRECT method for the inversion of the off-limb Stokes profiles indicated in the text (upper 
panels and bottom left panel). The number N of function evaluations is shown on the top 
of each panel. The presence of a very deep and shallow strip where the minimum is located 
makes it difficult to obtain it. The DIRECT method needs at least 1000 function evaluations 
to find the region where the minimum is located. The bottom right panel shows the value 
of the merit function with the white dot indicating the combination of parameters that 
give the smallest value of the merit function. This is an example of a problem that poses 
severe difficulties to any gradient-based method like Levenberg-Marquardt. 
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Fig. 19. — Value of 9b and h at which the merit function is evaluated when using the DIRECT 
method for the inversion of the on-disk Stokes profiles indicated in the text. The number 
N of function evaluations is shown on the top of each panel (upper panels and bottom left 
panel). Even with only = 50, the DIRECT method is able to locate the region of the 
global minimum. When the number of evaluations increases, it converges towards the values 
6b = 40° and h = 20". The bottom right panel shows the value of the merit function 
with the white dot indicating the combination of parameters that give the smallest value of 
the merit function. 
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Table 1: Parameters of the inversion. 



Free parameter Description Units 

B Magnetic field strength gauss 
6b Inclination of the magnetic field vector with respect to the vertical degrees 

Xb Azimuth of the magnetic field vector degrees 

Vth Thermal velocity affecting the width of the line km s~^ 

Vraac Bulk velocity of the plasma leading to a red/blue-shift km s~^ 

At Optical depth of the line — 

a Line damping parameter — 



Table 2: Inversion scheme. 



Step Method Free parameters Stokes profiles 

"1 DIRECT T;th, ^;mac, At, a I 

2 LM vth, fmac, At, a I 

3 DIRECT B, Ob, Xb I, Q, U, V 

4 LM B,9b,Xb I, Q, U, V 
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Table 3: He I atomic data. 



Transition {u — > £) Air wavelength [A] A{j3uLuSJu — > PiL^SJi) [s ^ 



[G] 



2p =^Po 


- 2s ^Si 


10829.0911 


1.022x10^ 


undefined 


2p 


- 2s ^Si 


10830.2501 


1.022x10^ 


0.77 


2p 


- 2s 


10830.3398 


1.022x10^ 


0.77 


3p =^Po 


- 2s ^Si 


3888.6046 


9.478x10^ 


undefined 


3p 


- 2s 


3888.6456 


9.478x10*^ 


0.72 


3p 


- 2s ^^Si 


3888.6489 


9.478x10^ 


0.72 


3s ^Si - 


- 2p 


7065.7085 


3.080x10^ 


0.18 


3s - 


- 2p 3Pi 


7065.2150 


9.250x10^ 


0.53 


3s ^Si - 


- 2p 3P2 


7065.1769 


1.540x10^ 


0.88 


3d ^Di 


- 2p ^Po 


5875.9663 


3.920x10^ 


8.92 


3d 3D2 


- 2p 


5875.6405 


5.290x10^ 


5.16 


3d 


- 2p 3Pi 


5875.6251 


2.940x10^ 


6.69 


3d ^Dg 


- 2p 3P2 


5875.6150 


7.060x10^ 


6.02 


3d 3D2 


- 2p 3P2 


5875.6141 


1.760x10^ 


1.72 


3d 


- 2p 3P2 


5875.5987 


1.960x10^ 


0.45 



